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Abstract 

A search for sleptons, neutralinos, charginos, sgoldstinos and heavy stable 
charged sleptons in the context of scenarios where the lightest supersymme- 
tric particle is the gravitino, is presented. Data collected during 2000 with the 
DELPHI detector at centre-of-mass energies from 204 to 208 GeV were analysed 
and combined with all the data collected from 1995 to 1999 at lower energies. 
No evidence for the production of sleptons, neutralinos and charginos has been 
found, therefore new limits on the mass of these supersymmetric particles and 
on the model parameter space are set. The search for heavy stable charged 
sleptons also updates the stable sleptons mass limit. The absence of evidence 
for sgoldstino production allows limits to be set on its mass and on the scale of 
supersymmetry breaking. 
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1 Introduction 

In the year 2000 the Large Electron Positron Collider (LEP) at CERN finished its ope- 
ration achieving record energies of 204 to 208 GeV when the DELPHI detector collected 
an integrated luminosity of 223.53 pb^^. These data were analysed to update the searches 
for sleptons, neutralinos, charginos, sgoldstinos and heavy stable charged particles [1,2] 
in the context of gauge mediated supersymmetry breaking (GMSB) models. 

Supersymmetry (SUSY) is usually assumed to be broken in a hidden sector of parti- 
cles and then communicated to the observable sector (where all the particles and their 
superpartners lie) via gravitational interactions. An alternative possibility is that this 
mediation is performed by Standard Model (SM) gauge interactions, leading to models of 
gauge mediated supersymmetry breaking. In most current GMSB theoretical work [3-5] , 
it is assumed that the hidden sector is coupled to a messenger sector, which in turn cou- 
ples to the visible sector through radiative corrections with gauge-interaction strength. 
The primary motivation for GMSB is that it naturally accommodates the experimentally 
observed absence of fiavour changing neutral currents due to the fact that gauge interac- 
tions arc fiavour blind. In these models the scale of SUSY breaking (\/F) can be as low 
as about 10"^ or 10^ GeV^ in order to have supersymmetric particle (sparticle) masses of 
the right order of magnitude (~ 100 GeV/c^). 

The mass of the gravitino (G) is related to the scale of SUSY breaking through the 
expression: 

ma ~ 2.5 x F/(100 TeV)^ eV/c^ (1) 

therefore can be as low as few eV/c^. Consequently in these models G is the lightest 
supersymmetric particle (LSP) and all the other sparticles will decay into final states 
that include it. Gravitino masses below 3-10"'^ eV/c^ have been ruled out using Tevatron 
data for the multijet final states [6]. On the other hand, hints from cosmology point in 
the direction of either a light gravitino with mass below 1 keV/c^ or a heavy one with 
mass above 1 TeV/c^ [7,10]. 

In GMSB models the entire minimal supersymmetric standard model (MSSM) spec- 
trum can be predicted in terms of the following parameters: 

F, A, M, n, tan/? and sign (/x). (2) 

The most important parameter is A (the effective SUSY breaking scale) because it sets 
the overall mass scale of supersymmetric particles. M is the messenger mass scale. The 
number of messenger generations, n, is also very important because it determines which 
sparticle is the next-to-lightest supersymmetric particle (NLSP). For n = 1 the NLSP is 
mainly the and for n > 2 it is one of the sleptons. The parameter tan/3 is the ratio of 
the Higgs vacuum expectation values, and sign(//) is the sign of the Higgs sector mixing 
parameter^. 

The coupling to the gravitino is very weak, therefore, all the superparticles other than 
the next-to-lightest supersymmetric particle undergo chain decay down to the NLSP 
which finally decays to the G. The mean decay length (L) of the NLSP depends on 
ma [8] . Namely, for I — > IG decay: 

L = 1.76 X lO-Jf^V-l ( -^^] " ( ' cm (3) 
V \mj VlOOGeV/cV \leV/cy ^ ' 

iln gravity mediated SUSY breaking models a/F ~ lO^" or lO" GeV. 

^The magnitude of jj, is calculable from the other parameters in the model by imposing radiative electroweaJj symmetry 
breaking. 
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where mf is the slepton mass and Ej the slepton energy. Therefore, the gravitino mass 
determines if the NLSP decays inside or outside the detector, giving rise to very interesting 
topologies explored in this paper. For example, for ttiq < 250 eV/c^{y/F < 1000 TeV), 
the decay of a NLSP with mass greater than for example 60 GeV/c^ can take place within 
the detector. This range of \/F is in fact consistent with astrophysical and cosmological 
considerations [9,10]. Figure 1 shows the f mean decay length as a function of the 
gravitino mass for different f masses. 

In this paper data were analysed within the two possible slepton NLSP scenarios as 
discussed in the following. Depending on the magnitude of the mixing in the third family 
between the left and right gauge cigcnstates, and f^, there arc two possible scenarios. 
If the mixing is large"\ fi (the lighter mass eigenstate) is the NLSP. However, if the 
mixing is negligible, fi is mainly right-handed [11] and almost mass degenerate with the 
other sleptons. In this case, the and flji three body decay (Z— > TitI with fi— > r G), is 
very suppressed, and cr and (Ir decay directly into IG. This scenario is called sleptons 
co-NLSP. 

The signature for SUSY particle production within GMSB models at LEP2 depends 
on the NLSP type and on its mean decay length, or equivalently, on the gravitino mass. 
The NLSP could be pair produced directly, or other sparticle production could lead to a 
cascade decay into the NLSP. The NLSP will decay into its non-SUSY partner and a G. 
Taking into account all these factors, the following topologies can be expected: 

• For tjIq below a few eV/c^, the NLSP decays in the vicinity of its production 
point, before the tracking devices of the detector, and different topologies can be 
expected. If the sleptons are pair produced the signature in the detector is the same 
as in the search for sleptons in gravity-mediated supersymmetry breaking scenarios 
(MSUGRA) with Xi LSP, i.e. two acoplanar^ leptons and missing energy [12] since 
the slepton decays into a lepton and a ncutralino that escapes detection. This pro- 
cess is topologically equivalent to the pair production of two sleptons and each of 
them decaying into a lepton and a gravitino (if uIq is below a few eV/c^). However, 
if neutralino pair production is kinematically allowed, the production cross-section 
can be larger than for / even if m^o > mj because of the f3'^ suppression factor of the 
scalar production cross-section. In this case the topology is given by four leptons 
and missing energy since each x5 decays into II, and the sleptons into IG. The mass 
region that can be inspected using this search can be complemented with the search 
for lightest neutralino pair production when the neutralino is the NLSP decaying 
into a photon and a gravitino [13]. 

• For ttIq between a few eV/c^ and a few hundred eV/c^ the NLSP has an intermediate 
mean decay length and it would decay in flight in some part of the detector volume. 
This creates well defined secondary vertices or kinks when the / is reconstructed by 
the tracking devices, or large impact parameter tracks if it is not. 

• For gravitino masses above few hundred eV/c^ the NLSP would be sufficiently long- 
lived to decay outside the detector giving rise to heavy stable charged particle sig- 
natures. 

In the GMSB parameter space where the Xi is the NLSP, the chargino is always much 
heavier than 100 GeV/c^ [14] and, therefore, cannot be produced at LEP. On the contrary, 
in the parameter space where the / is the NLSP there are regions where the xf is light 
enough to be produced [15]. In this case the topology for below a few eV/c^ is again 

"^In GMSB models large mixing occurs generally in regions of tan/3 > 10 or > 500 GeV/c^. 

■'Acoplanarity is defined as the complement of the angle between the projections of the two tracks onto the plane 
perpendicular to the beam. 
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two acoplanar Icptons and missing energy since each xf decays into lu and each / into IG. 
For uiq between a few eV/c^ and a few hundred eV/c^ the topologies are kinks or large 
impact parameter tracks and, for gravitino masses above a few hundred eV/c^, heavy 
stable charged particle signatures are expected. 

In this paper, the update of the search for heavy stable charged particles is also 
performed. This kind of particle is predicted not only in GMSB models but also in MSSM 
with a very small amount of R-parity violation [16], or with R-parity conservation if the 
mass difference between the LSP and the NLSP becomes very small (references in [17]). 
The typical signature of these events is two massive particles traversing the detector which 
do not produce Cherenkov radiation in DELPHI'S Ring Imaging CHcrcnkov (RICH) 
detectors, but high ionization losses in the Time Projection Chamber. Updated lower 
mass limits on heavy stable charged particles, under the assumption that the LSP is a 
charged slepton, are presented. 

Recently it has been pointed out [18] that an appropriate theory must also contain the 
supersymmetric partner of the goldstino, called the sgoldstino, which could be massive. 
In the minimal R-parity conserving model, as considered in [18], the effective theory at 
the weak scale contains two neutral scalar states: S, CP-even and P, CP-odd (from now 
on the two states will be labelled with the generic symbol since the discussion applies to 
both of them) . It must be pointed out that sgoldstinos have even R-parity, therefore they 
are not necessarily produced in pairs and their decay chains do not necessarily contain 
the LSP. The production of these supersymmetric particles may be relevant at LEP2 
energies in light gravitino scenarios. One of the most interesting production channels is 
the process e+e" 07 which depends on the mass (m^) and on The most relevant 
decay modes are — > 77 and — > gg. The corresponding branching ratios depend on 
the gaugino masses Mi, M2 and M3, and the total width is P ~ r(0 — * 77) -|-r(0 — >^ gg). 
In this paper two sets for these parameters are considered as suggested in [18]; they are 
listed in Table 1. The total width for a large interval of the parameter space is narrow 
(below a few GeV/c^), except for the region with small Vf where the production cross- 
section is expected to be very large. The two decay channels considered produce events 
with very different topologies. The channel — > 77 gives events with three high energetic 

S TTlP" 

photons, one of which has monochromatic energy {E^ — ^ large fraction of the 

parameter space where has a negligible width. Despite the lower decay branching 
ratio (4 and 11% for the two sets of Table 1, respectively), this final state is worth 
investigating because the main background source is the QED process e"'"e~ —>■ 77(7), 
which is expected to be small if photons in the forward region are discarded. On the 
other hand, the channel S ^ gg gives events with one monochromatic photon (except for 
the region with small \fF) and two jets. An irreducible background from e^e~ ggj 
events is associated to this topology and therefore the signal must be searched for as an 
excess of events over the background expectations for every mass hypothesis. 





Ml 


M2 


M3 


BR{(j) 77) 


BR{<P ^ gg) 


1) 


200 


300 


400 


4% 


96% 


2) 


350 


350 


350 


11% 


89% 



Table 1: Two choices for the gaugino mass parameters (in GeV/c^) relevant for the 
sgoldstino production and decay, and the corresponding branching ratios (BR) of the two 
channels considered. The BR are almost independent of m<^ in the mass region below 
200 GeV/c2. 
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The list of GMSB signatures analysed in this paper is given in Table 2. 



Production 


Decay mode 


L 


Expected topology 


e+e- ^ I! 


I^IG 


^detector 
^ ^detector 
^detector 


Acoplanar leptons 

Kinks and large impact parameters 

Heavy stable charged particles 


e+e- xSx? 


Xi^ll^ IIG 


^detector 


Four leptons 


e+e" ^ xixl 




^detector 

^detector 
^detector 


Acoplanar leptons 

Kinks and large impact parameters 

Heavy stable charged particles 


e+e" — > 07 


^ 77 


"^"^ ^detector 
■^"^ ^detector 


3 high energetic 7 

1 monochromatic 7 and 2 jets 



Table 2: Final state topologies studied in the different scenarios. 



The organization of the paper is as follows. A brief description of the DELPHI detector 
is presented in section 2. The data samples are described in section 3. The different 
selection criteria, the efficiencies and the number of events selected in data and in the 
expected Standard Model background are reported in section 4. Finally, the results are 
presented in section 5 comprising cross-section limits of the pair produced sparticles, 
lower mass limits and limits on the GMSB model parameters. 

2 Detector description 

DELPHI was one of the four detectors operating at the LEP collider from 1989 to 2000. 
It was designed as a general purpose detector for e+e" physics with special emphasis on 
precise tracking and vertex determination and on powerful particle identification. A 
detailed description of the DELPHI detector can be found in [19] and the detector and 
trigger performance in [20,21]. Here only those components relevant for the present 
analyses are discussed. 

Charged particle tracks are reconstructed by a system of tracking chambers inside the 
1.2 T solenoidal magnetic field: the Vertex Detector (VD), the Inner Detector (ID), the 
Time Projection Chamber (TPC) and the Outer Detector (OD) in the barrel region; two 
sets of planar drift chambers aligned perpendicular to the beam axis (Forward Chambers 
A and B) measure tracks in the forward and backward directions. 

For the data presented here, the VD consists of three cylindrical layers of sihcon 
detectors, at radii 6.3 cm, 9.0 cm and 11.0 cm, and polar angle acceptance from 24° to 
156°. All three layers measure coordinates in the plane transverse to the beam (xy), and 
at least two of the layers also measure z coordinates along the beam direction. The ID 
consists of a cylindrical drift chamber with inner radius 12 cm and outer radius 22 cm, 
surrounded by 5 layers of straw tubes, having a polar acceptance between 15° and 165°. 
The TPC, the principal tracking device of DELPHI, consists of a 2.7 m long cylinder of 
30 cm inner radius and 122 cm outer radius. Each end-plate of the TPC is divided into 
6 sectors with 192 sense wires and 16 circular pad rows per sector. The wires help in 
charged particle identification by measuring the specific energy loss (dE/dx) and the pad 
rows are used for 3 dimensional space-point reconstruction. The OD consists of 5 layers 
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of drift cells at radii between 192 cm and 208 cm, covering polar angles between 43° and 
137°. 

The electromagnetic calorimeters consist of a High Density Projection Chamber 
(HPC) covering the polar angle region from 40° to 140° and, a Forward ElcctroMag- 
netic Calorimeter (FEMC) covering the polar angle regions from 11° to 36° and from 
144° to 169°. The Scintillator Tile Calorimeter (STIC) extends the polar angle coverage 
down to 1.66° from the beam axis in both directions. The Hadron CALorimeter (HCAL) 
covers 98% of the solid angle. The muons which traverse the HCAL are recorded in a set 
of Muon Drift Chambers placed in the barrel, forward and backward regions. 

The Ring Imaging CHerenkov (RICH) detectors of DELPHI provide charged particle 
identification in both the barrel (BRICH) and forward (FRICH) regions. They contain 
two radiators of different refractive indices. The liquid radiator is used for particle iden- 
tification in the momentum range from 0.7 to 8 GeV/c. The gas radiator is used for 
particles with momentum range from 2.5 GeV/c to 25 GeV/c. 

3 Data sample and event generators 

The searches reported in this paper are based on data collected with the DELPHI 
detector during 2000 at centre-of-mass energies from around 204 to 208 GeV. The total 

integrated luminosity was 223.53 pb^^. Table 3 summarises the energies analysed and 
the integrated luminosities corresponding to each energy during the LEP2 period. 



Year 


1995 


1996 


1997 


1998 


1999 


2000 


(GeV) 


130-136 


161-172 


183 


189 


192-202 


204-208 


£ (pb-) 


11.9 


19.6 


54.0 


158.0 


228.2 


223.5 



Table 3: Centre-of-mass energies analysed and their corresponding integrated luminosities 
during the LEP2 period. 

To evaluate the signal efficiencies and background contamination, events were genera- 
ted using different programs, all relying on JETSET 7.4 [22], tuned to LEPl data [23] for 
quark fragmentation. 

Slepton pair samples at 208 GeV centre-of-mass energy were produced with PYTHIA 
5.7 [22] with sleptons having mean decay lengths from 0.25 to 200 cm and masses from 
60 to 104 GeV/c^. Other samples of slepton pairs were produced at 206 and 208 GeV 
with SUSYGEN [24] for the small impact parameter search with rrif from 90 GeV/c^ to 
102 GeV/c^ and equal to 90 GeV/c^. A sample of selectrons with mass equal to 
90 GeV/c^ was produced at 206 GeV. Neutralino pair events and their subsequent decay 
products were generated with SUSYGEN. Selection efficiencies were computed from samples 
with neutralino masses from 72 GcV/c^ < mj + 2 GeV/c^ < m^o < ^/s/2 at 206 GeV. 
SUSYGEN was also used to generate the chargino pair production and decay. In order to 
compute detection efficiencies, samples at 204 GeV and 206 GeV centre-of-mass energies 
were generated with gravitino masses of 1, 100 and 1000 eV/c^, m^-l-0.3 GeV/c^ < < 
\/s/2 and 80 GeV/c^ < mf < \/s/2 — 2.6 GeV/c^. Samples with smaller Am = m-+ —ruf^ 
were not generated because in that region each chargino decays into a W and a gravitino 
with an appreciable branching ratio. 

In the search for heavy stable charged particles, signal efficiencies were estimated from 
pair produced heavy smuons generated at energies of 205 GeV, 206.7 GeV and 208 GeV 
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with SUSYGEN. The events were passed through the detector simulation as heavy muons. 
The efficiencies were estimated for masses between 10 GeV/c^ and 100 GeV/c^. 

For the sgoldstino search, signal efficiencies for the channel — > 77 were estimated 
from QED background events generated according to [25]. On the other hand, the se- 
lection efficiency for the (j) ^ gg channel was evaluated using gg7 background events 
generated with PYTHIA and processed through the full DELPHI analysis chain and re- 
weighted according to the background and signal photon polar angle distribution. 

The background process e+c^ qq{nrf) was generated with PYTHIA 6 . 125, while 
KORALZ 4.2 [26] was used for /i+/i"(7) and r+r"(7). The generator BHWIDE [27] was 
used for e"'"e" — > e+e" events. Processes leading to four-fermion final states were gen- 
erated using EXCALIBUR 1.08 [28] and GRC4F [29]. Two-photon interactions leading to 
hadronic final states were generated using TWOGAM [30], including the VDM, QPM and 
QCD components. The generators of Berends, Daverveldt and Kleiss [31] were used for 
the leptonic final states. 

The cosmic radiation background was studied using cosmic muons collected before the 
beginning of the 2000 LEP run. 

The generated signal and background events were passed through the detailed simu- 
lation [20] of the DELPHI detector and then processed with the same reconstruction and 
analysis programs used for real data. 

4 Data selection 

The following sections describe the selection criteria used to search for the different 
topologies summarized in Table 2. The searches for slepton, neutralino and chargino pair 
production are based on the ones described in [1,32-35]. The search for heavy stable 
charged sleptons is based on the analysis already published in [1,36,37]. Finally the 
sgoldstino search has already been presented in [2] . 

4.1 Slepton pair production 

This section describes the selection criteria used in the search for the process e+e~ —>■ 
I'^l" l^Gl~G. Loose preselection cuts were imposed on the events in order to suppress 
as much as possible the low energy background (beam-gas and beam-wall) and the SM 
processes. The reconstructed tracks of charged particles were required to satisfy certain 
quality criteria: momenta above 100 MeV/c and impact parameters (geometric signed) 
below 4 cm in the plane transverse to the beam pipe (xy), and below 10 cm in the direction 
along the beam pipe {z). Clusters in the calorimeters were interpreted as neutral particles 
if they were not associated to charged particles and if their energy exceeded 100 MeV. 
Only these particles were used to compute the general event quantities. The preselection 
cuts were the following: 

• to eliminate high multiplicity events like e+e^ (HijT'l) 01 WW, Wez/g and ZZ 
when the produced particles had pure hadronic or semileptonic decays, the multiplic- 
ity computed with the particles that satisfied the quality requirements was required 
to be between 1 and 6 (the multiplicity of all signal samples was very well contained 
between these two limits); 

• to eliminate two-photon processes, the visible energy in the event was required to 
be above 10 GeV; 
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• to eliminate the remaining contribution of two-photon and two-fermion processes, 
the absolute value of the transverse momentum vector of charged and neutral par- 
ticles was required to be greater than 5 GeV/c; 

• the energy measured in the very forward calorimeters (STIC) was required to be 
below 10 GeV to eliminate the residual contamination of processes mentioned above; 

• to eliminate the Bhabha contribution, the total electromagnetic energy was required 
to be less than the beam energy. 

All the events that survive the preselection cuts underwent the search for secondary 
vertices or kinks. Only the events which were not tagged as kink candidates passed the 
selection criteria to search for large impact parameter tracks. 

4.1.1 Secirch for secondciry vertices or kinks 

The analysis exploits a peculiarity of the l^G topology in the case of in- 

termediate gravitino masses (i.e. few eV/c^< niQ < few hundred eV/c^), namely, one 
or two tracks coming from the interaction point and at least one of them with either a 
secondary vertex or a kink. 

All the charged particles of the event that survived the preselection cuts were grouped 
into clusters (in order to group all the particles coming from a tau decay) according to 
their first measured point in the xy plane. This clustering procedure was iterative and 
worked as follows. The pair of particles with the smallest separation at their respective 
starting points was considered first. If this separation was smaller than 2 cm, the par- 
ticles were grouped to form a cluster whose starting point was defined as the average of 
their first measured points. The two particles were then replaced by this cluster which 
was subsequently treated as a pseudo-particle. The process was then repeated until all 
charged particles or pseudo-particles were grouped into clusters. This procedure allowed 
for clusters containing a single particle if its momentum was larger than 1.5 GeV/c. 
Events were rejected if more than 6 particles were not grouped into clusters or if a cluster 
could not be obtained. This cut was intended to eliminate the remaining beam related 
background events^ that had not been excluded at the preselection level. 

Once all the particles were grouped in clusters, the search for kinks was performed 
in the following way. Slepton candidates were searched for among all the clusters in the 
event. Among the remaining clusters, the ones corresponding to lepton candidates or 
decay products of taus were also searched for. The clusters were extrapolated in order 
to find a crossing point. If the crossing point existed, the event was considered as a kink 
candidate. Reconstruction of secondary vertices for the case f — > tG is illustrated in 
Figure 2, which shows a decay vertex and the variables used in the analysis.^ 

Isolated particles (clusters with only one particle) were considered as / candidates 
if their trajectories were compatible with particles coming from the interaction point 
according to the following selection criteria: 

• the first measured point with respect to the beam spot in the plane transverse to 
the beam axis {R^p) had to lie in the Vertex Detector (VD); 

• the momentum of the particle was greater than 2 GeV/c; 

• the polar angle with respect to the beam axis had to satisfy |cos^| < 0.8, corres- 
ponding to the barrel region; 

• the impact parameter along the beam axis and in the plane perpendicular to it was 
less than 10 and 4 cm, respectively. 

^This kind of event was mainly characterized by a high number of low momentum charged particles seen only in the 
innermost detectors VD and ID. 
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For every / candidate, a search was made for a second cluster satisfying the following 
selection criteria: 

• the starting point in the transverse plane (-R^p) ^'^ greater than i?^^. The 
second cluster starting point was always found in the Inner Detector (ID) or the 
Time Projection Chamber (TPC); 

• the angular separation between the directions defined by the I candidate and the 
lepton candidate had to be smaller than 45° in the xy plane, to consider only the 
particles which were in the direction of the slepton. 

The / candidate and the lepton cluster had to define a common crossing point, called 
secondary vertex or kink. If the lepton cluster included more than one charged particle 
(which is the case when Hs a r decaying to 3 or 5 prongs) , only the one with the highest 
momentum was used to search for the kink. To find the crossing point, the particle tra- 
jectories were represented by a helix in space. Taking into account this parametrization, 
the point of closest approach between the / particle and the selected particle from the 
Id cluster was calculated. The conditions to define a good crossing point between both 
particles were the following: 

• the minimum distance between the particles had to be smaller than 1 mm in the xy 
plane; 

• the crossing point, the end point of the slepton and the starting point of the lepton 
were required to satisfy the following conditions: 

-10 cm < {Rcross - Rind) < 25 Cm 

-25 cm < (Rcross-Rip) < 10 cm, (4) 

where R^nd' Rcross and are the distances w.r.t the beam spot of the end point of 
the slepton, the crossing point and the starting point of the lepton in the xy plane. 
The cut was optimised to assure that all VD or IDVD only particles (particles which 
only had ID or VD hits) were connected with particles reconstructed with the TPC. 

The resolution achieved (generated distance minus reconstructed distance) with the algo- 
rithm to find secondary vertices in the coordinates x and y was 0.14 mm. The z coordinate 
was not taken into account in the search for a crossing point because the resolution for 
low momentum particles is very poor. 

Fake decay vertices could be present among the reconstructed secondary vertices. They 
could be produced by particles interacting in the detector material or by radiated photons 
when the particle trajectory was reconstructed as two separated particles. To eliminate 
these events, additional conditions were required: 

• to reject hadronic interactions, the angle between the direction of any reconstructed 
hadronic vertex w.r.t. beam spot (secondary vertices reconstructed in region where 
there is material) and the direction of the slepton candidate must be greater than 

5°; 

• to reject segmented tracks, the angle between the particles used to define a vertex 
had to be larger than 6°; 

• to reject photon radiation in the case of / clusters with only one particle, there had 
to be no neutral particle in a 3° cone around the direction defined by the difference 
between the I momentum and the momentum of the I calculated at the crossing 
point. 
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If no pair of particles was found to survive these conditions, the event was rejected. 
Figure 3 shows the distribution of these three quantities. The distributions compare real 
data, expected SM background simulation and a simulated signal for nif— 60 GeV/c^ with 
a mean decay length of 50 cm. 

Table 4 shows the different SM background contributions and the observed events in 
data after applying the selection criteria to search for kinks. Efficiencies for different 
gravitino and stau masses were calculated by applying the above selections to the sim- 
ulated signal samples. Figure 4 shows the secondary vertex reconstruction efficiency as 
a function of the stau decay radius. For smuons and selectrons the same dependence 
is observed. For smuons the efficiency plateau is around 60%, while for selectrons it is 
around 40% due to the preselection cut on total electromagnetic energy. 



Observed events 


2 


Total background 


o.88^i:- 




lfi+*J-U9 
U--LU_o.05 




00+'^-*'' 


4-fermion (except 77) 


12+"-"'^ 
"-"•-"-^-o.oi 


77 — >• T^T~' 


'^•-'-"-0.06 


77 — > e^e~ 


0.44l[j;^^« 



Table 4: Number of observed events at i/i from 204 GeV to 208 GeV together with 
the total number of expected SM background events and the expected numbers from the 
individual background sources, for the secondary vertex search. The asymmetric errors 
are due to the Poissonian description of the statistics. 

4.1.2 Lcirge impact pEirameter seeirch 

To investigate the region of low gravitino masses (short decay lengths) the previous 

search was extended to the case of sleptons decaying between 0.25 cm and around 10 cm, 
i.e., before the tracking devices. In this case it was only possible to reconstruct the slcpton 
decay products. The impact parameter search was only applied to those events accepted 
by the same preselection cuts as in the search for secondary vertices but not selected by 
the vertex analysis. The events used in this search contained exactly two single particle 
clusters (i.e. two charged particles with momentum larger than 1.5 GeV/c and a distance 
between starting points greater than 2 cm) which were acoUinear and had large impact 
parameters. The events were accepted as candidates if: 

• the first measured point in the xy plane of at least one of the particles was in the 
VD; 

• both particles were reconstructed with the TPC to guarantee good particle recons- 
truction quality; 

• at least one of the particles had an impact parameter w.r.t. the beam spot larger 
than 0.2 cm in the xy plane to remove SM events; 

• the ratio of the maximum impact parameter over the minimum impact parameter 
in the xy plane and w.r.t. the beam spot was smaller than -1.5 or larger than -0.5, 
to reject cosmic rays since they are characterized by large impact parameters of the 
same value and opposite sign. 
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• the acoUinearity between the two particles was larger than 10° to eliminate back- 
to-back events with badly reconstructed particles or interactions which always gave 
small acoUinearities. In addition, to reduce further the cosmic ray background, 
the acoUinearity between the two particles was required to be smaller than 175°, 
since an off-time cosmic ray crossing from one TPC drift half to the other could be 
reconstructed as two almost parallel particles. 

Figure 5 shows the acoUinearity distribution for real data minus simulated SM back- 
ground, compared to cosmic rays for events which passed the preselection cuts and cuts 
1-3. The latter follows pretty well the acoUinearity distribution of the former difference. 

Tabic 5 shows the different SM background contributions and the number of events 
observed in data after applying the selection criteria to search for large impact parameter 
particles. 



Observed events 


2 


Total background 






'J-'J'^-o.oi 


Z77^(ee)(n7) 


n 1 9+o.yu 

U--L^_0.10 


4-fermion (except 77) 




77 — >• T^T~ 


u-uy_o.o4 


77 — >• e^e" 


1 4q+i-U8 

-■-•^^-0.34 



Table 5: Number of observed events at \fs from 204 GeV to 208 GeV together with 
the total number of expected SM background events and the expected numbers from 
the individual background sources, for the large impact parameter track search. The 
asymmetric errors are due to the Poissonian description of the statistics. 

The efficiencies were derived for different slepton masses and decay lengths by applying 
the same selection criteria to the simulated signal events. In the search for fi the maxi- 
mum efficiency was around 32% corresponding to a mean decay length of 2.5 cm. The 
efficiency decreased sharply for lower decay lengths due to the requirement on mini- 
mum impact parameter. For longer decay lengths, the appearance of reconstructed I 
in combination with the cut on the maximum number of charged particles in the event 
caused the efficiency to decrease smoothly. This decrease was compensated by a rising 
efficiency in the search for secondary vertices. For masses above 60 GeV/c^ no dependence 
on the / mass was found far from the kinematic limit. 

The same selection was applied to smuons and selectrons. For smuons the efficiency 
increased to ~ 58% for a mean decay length of 2.5 cm and masses over 60 GeV/c^ since 
the smuon always had a one-prong decay. For selectrons the efficiency was ~ 33% for the 
same mean decay length and range of masses. 

4.1.3 Small impact peirameter secirch 

The large impact parameter search can be extended further to mean decay lengths 
below 0.1 cm. Charged particles were selected if their impact parameter was less than 
10 cm in the plane transverse to the beam direction and less than 15 cm in the direction 
along the beam pipe. The polar angle had to be between 20° and 160°. Their measured 

''The a<;ollinearity is defined as being 180° minus the angle between the momentum vectors of both particles. 
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momentum was required to be larger than 400 MeV/c with relative error less than 100% 
and track length larger than 30 cm. Any calorimetric deposit associated to a discarded 
charged particle was assumed to come from a neutral particle. 

The search was restricted to events with 2 to 4 charged particles and missing energy 
larger than O.Sy^. The 77 events were suppressed by requiring a visible energy greater 
than 0.08^/s and a transverse missing momentum greater than 0.03\/i. The polar angle 
of the missing momentum was required to be between 30° and 150°, and the total energy 
in a cone of 30°, £"30, around the beam- pipe was required to be less than 10% of the 
visible energy, and the neutral energy was required to be less than 0.175y^. 

The events were then divided into two hemispheres using the thrust axis. The total 
momentum of charged and neutral particles in each hemisphere was computed and used to 
define the acoUinearity of the event. Standard e+e~ — > ffi'j) processes and cosmic rays 
were reduced by requiring the acoUinearity to be greater than 10°. The charged particle 
with the largest and good quality momentum {Api/pi < 50%) in each hemisphere was 
selected as the leading particle. The following quality requirements were only applied 
to the leading particles: the first measured point of the particle tracks had to be within 
50 cm of the beam spot in the xy plane, the particles were required to have at least 
one segment beyond the ID detector and to be away from insensitive regions of the 
electromagnetic calorimeter. In addition, at least one of the leading particles was required 
to be reconstructed with the TPC. 

e'^e~ — > ffi'j) processes and cosmic rays were further reduced by requiring an angle 
between the leading particles in the xy plane of less than 3 radians. Hadronic events, 
in particular 77 — > qq or, in general, events where the available energy is shared by 
many particles (including undetected/unselected), were rejected by requiring A/pf+~p| > 
0.03i/i, where pi and p2 are the momenta of the leading particles. To reduce Bhabha 
events the total electromagnetic energy of the leading particles, Ei + E2, had to be less 
than 0.35y^. By requiring that any leading particle with an impact parameter larger 
than 1 cm in the xy plane and measured w.r.t. the beam spot should be reconstructed 
by the TPC and at least one other detector, the residual cosmic rays (in particular the 
out-of-time cosmic rays) were rejected. Finally, photon conversion events with only two 
particles were rejected by requiring the angle between them at their perigee to be greater 
than 5°. 

The background left after the selection described above consisted mainly of events 
containing r pairs in the final state (7*/^* tt and WW — > tutu). To reject these 
events, the variable ^/bf + bl, where bi and b2 are the impact parameters in xy (measured 
w.r.t. the beam spot) of the two leading particles, was used. Requiring -y/fef+Tf > 
0.06 cm eliminated most of the remaining background. 

In order to preserve the efficiency in the region of decay length above 10 cm, where the 
I can be observed as a particle coming from the primary vertex and badly measured due 
to its limited length, further requirements on the particle quality were applied only to 
the leading particle with the largest impact parameter (measured w.r.t. the beam spot). 
This particle was required to have a relative momentum error < 30% and the particle had 
to be measured at least either in the TPC or in all of the other three tracking detectors 
in the barrel (VD, ID and OD). 

The efficiency of the search did not show any significant dependence on the / mass for 
masses over 40 GcV/c^ and far from the kinematic limit, and it could be parameterized as 
a function of the / decay length in the laboratory system. The efficiency for fi detection 
reaches ~ 40% for decay lengths around 2 cm. It is still 16% for a decay length of 0.1 cm. 
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and 13% for 20 cm. The efficiency for /j, detection reaches 45% around 2 cm, 15% at 
0.1 cm, and 17% at 20 cm. 

In order to increase the efficiency in the search for selectrons, the cut 
(£^1 + E2) < 0.35^/s was not apphed. The Bhabha events that survived the selection 
were those where at least one of the electrons underwent a secondary interaction, thus 
acquiring a large impact parameter. However, it was found that in these cases the mea- 
sured momentum of the electron was smaller than the electromagnetic energy deposition 
around the electron track. Therefore, the cut {Ei/pi + E2/P2) < 2.2 was used for the 
selectron search. The maximum efficiency reached in the selectron search was ~ 35% at 
~ 2 cm mean decay length. 

The number of events selected in the data was 4 in the f and pt, search. The same 
4 events also passed the e search. The expected SM background in both searches was 
3.3±0.3 events. Figure 6 shows the \/h\ + 6| distribution for data (dots), simulated 
SM backgrounds (grey histogram) and simulated signal of irif^ = 90 GeV/c^ and rriQ = 
25 eV/c^ nt ^/s = 206 GeV and a boosted mean decay length of around 1 cm (white 
histogram) after all other cuts. This figure only shows two of the candidates. The other 
two events not shown have v^&f-T&f ~ 1.5 cm. The overflow bin has 1.3 expected events 
from Monte Carlo and 2 from data. All the selected candidates were compatible with SM 
events. 

4.2 Neutralino pair production 

In this section, the selections used to search for the process e~^e~ — > XiXi ~^ 
fiTTiT — >• tGttGt within the fi NLSP scenario, and the process e+e^ — >• XiXi ~^ 
lll'V IGll'GV (with BR(x? ^ //) = 1/3 for each leptonic flavour) within the co-NLSP 
scenario, are presented. 

In the following, the preselection of the events, common to both scenarios, is pre- 
sented. The reconstructed tracks of charged particles were required to have momenta 
above 100 MeV/c and impact parameters below 4 cm in the transverse plane and be- 
low 10 cm in the longitudinal direction. The relative error on the measurement of the 
momentum had to be smaller than 100%. Clusters in the calorimeters were interpreted 
as neutral particles if they were not associated to charged particles and if their energy 
exceeded 100 MeV. All charged and neutral particles that satisfy these criteria were con- 
sidered good particles and they were used to compute the relevant event quantities. To 
assure good quality of the data, the ratio of good to total number of particles was re- 
quired to be above 0.7. Particles that did not pass quality selection but had an associated 
calorimetric energy of at least 2 GeV had their angles taken from those of the track, but 
their momentum was recomputed from the energy of the calorimetric measurement (such 
particles were not included in the good sample). Events had to have between four and 
ten good charged particle tracks. In addition, it was required that the thrust be less than 
0.99; the transverse momentum had to be bigger than 3 GeV/c, and | cos^^p^^^^ J < 0.95 
(polar angle of the missing momentum vector). Very forward-going events were elimi- 
nated by requiring E30 less than 70% of the total visible energy. With this preselection, 
the total number of simulated background events and real data events was reduced by a 
factor of about 6000. Only events passing these pre-selections were analysed further. 

The selection takes advantage of the fact that signal events can be separated into 
two different kinematic regions of the ijn^o^ ,mj) space: when the mass difference Am — 
m^o — nil bigger than about 10 GeV/c^, all four leptons carry similar momenta. When 
the difference becomes smaller, the two leptons coming from the decay of the I tend to be 
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the most energetic, increasingly so as the Xi mass increases. The Durham algorithm [38] 
was used to divide the event into four jets by allowing the jet resolution parameter to 
vary as a free variable. Numbering the jets from 1 to 4 with Ei > E2 > E3 > E4, a 
variable r was defined as: 

E3 + E4 , . 

r = — — - . (5) 

El + E2 

The distribution of r shifts towards lower values with increasing neutralino masses. 

At the preselection level the two main differences between the fi NLSP and Ir co- 
NLSP scenarios come from the fact that the mean number of neutrinos carrying away 
undetected energy and momentum and the number of charged particles per event is 
considerably bigger for the former scenario. 

In the fi NLSP scenario, the simulated background samples were then divided into 
two samples above and below r = 0.1 and different requirements were imposed in the 
two cases. No significant dependence on this variable was observed in the Ir co-NLSP 
scenario. Two sets of cuts were applied in order to reduce the 77 and ff (7) backgrounds 
and a third set of cuts to select events according to their topology. Those cuts are 
compiled in Table 6. 







fi NLSP 


Ir co-NLSP 






r > 0.1 


r < 0.1 




Cuts 


Et (GeV) 


> 11 


> 12 


> 4 


against 


^30 (GeV) 




< 


< 0.6y/s 


77 


rrimiss (GeV) 


< 0.88y/s 


< 0.9y/s 


< 0.88y/s 




Pmax ( GeV/c) 


> 4 


> 3 


> 8 




Pt ( GeV/c) 






> 6 


Cuts 


■^good 


< 7 


< 9 


< 7 


against 


T 


< 0.975 


< 0.975 


< 0.95 


ff(7) 


acoplanarity (°) 


> 8 


> 8 


> 8 




m„,iss (GeV) 


> 0.3v^ 




> 0.2v/i 


Cuts 


jet-beam angle (°) 


> 17 


> 17 


> 18 


based on 


2-jets cone angle (°) 


> 20 


> 20 


> 25 


topology 


4-'J<''t''^ scpaiutimi (°) 


> 8 


> 1 


> 9 



Table 6: Sets of cuts applied in the search for neutralino pair production to eliminate the 
different background sources. Pmax is the momentum of the charged particle with largest 
momentum; Ngood is the number of good particles in the event; T is the event thrust; 
rrirniss is the missing mass of the event. Signal events tend naturally to cluster into a 
4-jet topology. Taking this into account, cuts based on topology were applied: jet-beam 
angle is the angle between the jet direction and the beam direction; 2-jets cone angle 
represents the broadness of the 2-jets (considering only charged particles) obtained when 
reducing the 4-jct topology into a 2-jets configuration with the Durham algorithm; and 
finally, .^-jete separation is the angle between jets. 



After these cuts, an efficiency between 26 and 44% was obtained for the signal events 
in the fi NLSP scenario, and between 35 and 46% in the / co-NLSP scenario. The number 
of events remaining in data and simulated samples after the selection procedure were 8 
and 7.1±0.6 respectively in the fi NLSP scenario, and 7 and 6.6±0.6 respectively in the 
I co-NLSP scenario. 
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4.3 Chargino pair production 

The search for chargino pair production, e"'"e~ XiXi ~^ l^vl^V l^uGl^vG, 
makes use, without modification, of four different analyses depending on the gravitino 
mass or, equivalently, on the mean decay length of the slepton. When the slepton decays 
at the vertex, the combination of two analyses can be exploited, the search for charginos 
and the search for acoplanar leptons in gravity mediated supersymmctry breaking sce- 
narios. Details of these analyses can be found in [39,40]. For intermediate mean decay 
lengths of the slepton the topology is large impact parameter tracks or kinks, therefore, 
the two analyses explained in sections 4.1.1 and 4.1.2 can be used. Finally, if the slepton 
decays outside the tracking devices the signature corresponds to stable heavy leptons and 
this analysis is explained in section 4.4. 

The selections developed for these searches were thus applied to simulated data samples 
with different gravitino masses, and the results are presented in terms of 95% confidence 
level (CL) excluded regions in the (m^-, m-+) plane in section 5.3. 

4.4 Heavy stable charged particles search 

In this analysis it is assumed that the slepton lifetime is large enough that the 
sleptons can pass through the tracking devices without decaying. If the slcptons are 
pair produced, one expects two slowly moving particles, which are characterized by an 
abnormally large energy loss in the TPC Furthermore, no Cherenkov light is expected 
to be produced by them in the Ring Imaging CHerenkov detectors (RICH) of DELPHI. 
Both set of information are combined to identify such particles. 

Only events with two or three charged particles were considered. Candidate particles 
were required to have hits both in the VD and in the TPC. The particle momentum had to 
be above 5 GeV/c^, and the particle length was required to be at least 30 cm. To reduce 
cosmic rays background the absolute value of the particle impact parameter was required 
to be below 0.15 cm in the xy plane and 1.5 cm in z. Furthermore, the difference of the 
impact parameters in z of the two most energetic particles in the event was required to be 
less than 1 cm. An additional protection was added against showering electrons for which 
the particle extrapolation through the RICH can be unrehable. For charged particles with 
an associated electromagnetic energy exceeding 15 GeV, an additional hit in the Outer 
detector was required. A careful run selection ensured that the RICH detectors were fully 
operational. 

The energy loss measurement was required to be based on at least 80 wires, and the 
dE/dx cahbration was checked on a TPC sector by sector basis, using Z° calibration data 
of the same year. 

Events were selected if they contained at least one charged particle with: 

(I) momentum above 5 GeV/c, high ionization loss in the TPC and no associated 7s in 
the gas radiator of the RICH (gas veto) or, 

(II) momentum above 15 GeV/c, ionization loss at least 0.3 below the expectation for a 
proton and surviving the gas veto or, 

(III) momentum above 15 GeV/c, surviving the gas and the liquid RICH veto. 

An event was also selected if the two particles with the highest momentum were charac- 
terized by a high ionization loss or a gas veto, or if both particles had a low ionization 
loss. Special care has been taken about the dE/dx in sector 6 (S6) of the Time Projec- 
tion Chamber which was not operational during the second half of data taking in 2000, 
reducing the efficiency of the search by several percent as the dE/dx search windows 
could not be applied to particles pointing to S6. To recover some sensitivity for particles 
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pointing to this sector, only the double veto search window was applied, requiring hits in 
the Vertex Detector, Inner Detector and Outer Detector to ensure a good propagation of 
the particles through the RICH. 

The data of the year 2000 has been subdivided into 3 energy bins, corresponding to 
energies below 206 GeV (85 pb^^), between 206 GeV and 207 GeV (124 pb^^), and above 
207 GeV (11.4 pb~^). A total background of 0.25±0.04 events was estimated from data 
itself by counting the number of particles passing the individual selection criteria. 

No candidate events were selected in data. Figure 7 shows the data and the three 
main search windows. The expectation for a 95 GeV/c^ mass signal is also shown. Signal 
efficiencies were estimated from simulation. For particle masses between 10 GeV/c^ and 
60 GeV/c^ the signal efficiencies are of the order of 30%. For larger masses they rise with 
increasing mass to about 76-78%. Then the efficiency drops rapidly when approaching 
the kinematic limit, and it is assumed to be zero at the kinematic limit. 

4.5 Sgoldstino search 

This section describes the search for e^e~ — > 07 events, with the sgoldstino 
going to two gammas or two gluons. The two channels considered here give rise to two 
different topologies. On the one hand, if the sgoldstino decays into two photons, the final 
topology of the event is three high energy photons, one of them monochromatic. On the 
other hand, if the sgoldstino decays into two gluons, one monochromatic photon and two 
jets can be expected in the final state. 

4.5.1 ^ 77 channel 

Events were selected as 777 candidates if they had at least two electromagnetic 
energy clusters with 0.219 < E/^ < 0.713; at least another one with E > 5 GeV and 
no more than two additional clusters, the second (if present) with E < 5 GeV. The two 
most energetic electromagnetic clusters had to be in the HPC region, 42° < 9 < 89° 
(91° <e < 138°), or in the FEMC region, 25° < ^ < 32.4° (147.6° <e < 155°). Finally, 
the third cluster had to be in the region 9 > 42° {9 < 138°) or 20° < ^ < 35° (145° < 
9 < 160°). The event should not have hits in two of the three Vertex Detector layers 
compatible (within ±2° in the azimuthal direction) with the extrapolated trajectory of a 
particle from the beam crossing point to an electromagnetic cluster in the calorimeters. 

Further, two hemispheres were defined by a plane orthogonal to the direction of the 
most energetic cluster. One hemisphere was required to have no charged particles de- 
tected in the barrel region of the tracking devices other than the VD with a momentum 
greater than 1 GeV/c extrapolating to within 5 cm of the mean beam crossing point. 
The requirement was further strengthened to suppress the larger e+e" background, by 
demanding that both hemispheres had no such particle detected by the TPC with 9 < 35°. 

The events obtained after this selection had a three-body final state kinematics if 
there was no significant initial state radiation lost along the beam pipe. Defining A = 
l^ul + l^isl + I ^23 1 1 where Sij is the angle between the particles i and j, A should be 
360° in a three-body final state since the particles lie in a plane. If only the events with 
A > 358° were accepted, the energies of the particles could then be determined with very 
good precision on the basis of the measured photon directions: 
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with S = sin 5i2+sin 5i3+sin 523. The error on the energy evaluation was further minimised 
by requiring min(5i2, ^13, ^23) > 2°. 

In 07 events the decay products are expected to be isotropically distributed in the 
(f) centre- of- mass system. This fact implies that the distribution of cosa, where a is the 
angle between the (f) direction (opposite to the prompt photon) and the direction of one 
of the two 4> decay products, in the centre-of-mass system, should be fiat. On the other 
hand, in the QED background, |cosa| peaks at 1 and therefore only the combinations 
giving Icosal < 0.9 were accepted. 

The number of selected events giving up to three combinations and the expected 
background were 22 and 20.31^:9, respectively. The error on the background was due to 
a correction applied in order to take into account the missing higher orders (additional 
radiation above which gives events having low values of A) in the simulation of the 
QED background. These kind of events were removed only from the selected sample 
of real data and therefore, a corresponding normalization correction factor of (—13^7)% 
was applied to the simulated sample. This correction is the dominant contribution to 
the systematic error. No significant background in addition to e"^e~ 77(7) events was 
found. 

No significant variation in the acceptance for a 07 signal and in the selection efficiency 
inside the acceptance region were observed in the 2000 data with respect to the lower 
energies in previous years: acceptance (51 ± 2)% and efficiency (76.6 ± 2.5)%. 

The energy resolution remained also unchanged and it was better than 0.5% over the 
whole photon energy range. 

The photon recoil mass spectrum obtained for the events collected during the 2000 
run including those reported in [2] is shown in Figure 8-a. The data are superimposed 
on the expected QED background distribution. 

4.5.2 (f) ^ gg channel 

This channel is expected to give a final state with one photon and two jets. An event 

was selected as a "-ygg candidate if it had an electromagnetic energy cluster identified 
as photon with > 5 GeV and 9 > 20°. The event must not have electromagnetic 
clusters below 9 = 5°. The total multiplicity (charged and neutral) had to be greater 
than 10, and the charged multiplicity greater than 5. To remove 77 events the cut 

Y17=i '\J (pI '^Py)i > 0.125-v/i (where n is the total multiphcity) was used. The sum of 
the absolute values of all particle momenta along the thrust axis had to be greater than 
0.20v^. An electromagnetic cluster with E < OAby/s, or a total particle multiplicity 
greater than 16 when the cluster energy was greater than 0.45y^ had to be present. 
The polar angle of the missing momentum had to satisfy \cos{9p^^^J\ < 0.995. The 
visible energy had to be greater than 0.60i/i. The jets had to be incompatible with 
the bb hypothesis by requiring the combined btag of the events to be less than zero [41]. 
Finally, the aforementioned | cosq;| and A had to be less than 0.9 and greater than 350°, 
respectively. 

The events were reconstructed forcing all particles but the photon into a 2-jet topology 
using the Durham [38] algorithm. Events were removed if Ucut > 0.02. The events were 
also rejected if the angle between the photon and the nearest jet was less than 10°. In 
the case of more than one photon candidate in the event, the most energetic one was 
considered as the one produced in e'^e~ —>■ 07. 

Similar to the 777 selection, the events obtained after this selection were three-body 
final state events in absence of additional lost radiation. Therefore kinematic constraints 
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were applied here as well. In this case, however, the jet direction was determined with a 
poorer precision than that obtained for photons, therefore the cut in A was less stringent 
and the resolution for the reconstructed photon energy was poorer: a two-Gaussian fit 
gave (Ji = 1.2 GeV (55% of the area) and (T2 = 4.1 GcV. 

The number of selected events and the expected background were 766 and 775 ± 5, 
respectively. 

No significant variation in the acceptance and in the selection efficiency inside the 
acceptance region was observed in the 2000 data compared to the values at lower energies. 
The acceptance was (76 ±2)% (almost independent of m,^) and the efficiency ranged from 
20 to 55% depending on the photon energy. The energy resolution was also unchanged. 

The photon recoil mass spectrum obtained for the events collected during the 2000 run 
and including those reported in [2] is shown in Figure 8-b. The data are superimposed 
on the expected background distribution. 

5 Results and interpretation 

Since there was no evidence for a signal above the expected background, the number 
of candidates in data and the expected number of background events were used to set 
hmits at the 95% confidence level (CL) on the pair production cross-section and masses 
of the sparticles searched for. The model described in reference [42] was used to derive 
limits within the GMSB scenarios. This model assumes radiatively broken clcctrowcak 
symmetry and null trilinear couplings at the messenger scale. The corresponding para- 
meter space was scanned as follows: 1 < n < 4, 5 TeV < A < 90 TeV, 1.1 < M/A < 10^ 
1.1 < tan/9 < 50, and sign^fi) = ± 1. The meaning of the parameters is explained in 
section 1. The limits presented here are at ^/s = 208 GeV after combining these results 
with those of the searches at lower centre-of-mass energies using the likelihood ratio 
method [43]. 

5.1 Slepton pair production 

The results of the search for slepton pair production are presented in the {mQ^mi) 
plane in Figure 9-a combining the impact parameter, the kink and the stable heavy 
lepton analyses, and using all DELPHI data from 130 GeV to 208 GeV centre-of-mass 

energies [1,32-34]. 

The fi pair production cross-section depends on the mixing in the stau sector. There- 
fore, in order to put limits on the fi mass, the mixing angle had to be fixed. The results 
presented here correspond to a mixing angle in the stau sector which gives the minimum 
fi pair production cross-section. Within the fi NLSP scenario, the impact parameter 
and kink analyses extended the limit m^^ > 82.5 GeV/c^ for itlq < 6 eV/c^, set by 
MSUGRA searches [12], up to tuq = 600 eV/c^, reaching the maximum excluded value 
of rufi— 93.6 GeV/c^ for ttiq— 130 eV/c^. For uiq > 130 eV/c^ the best lower mass hmit 
was set by the stable heavy lepton search. Although not shown in the figures, DELPHI 
excludes rriQ from 3-10~^ eV/c^ to 0.2 eV/c^. nif^ masses below 40 GeV/c^ have been 
already excluded as discussed in [32]. 

Within the sleptons co-NLSP scenario, the cross-section limits were used to derive 
lower limits for (Figure 9-b) masses at 95% GL. Assuming mass degeneracy between the 
sleptons, the kink and impact parameter searches extended the limit > 88 GeV/c^ set 
by MSUGRA searches [12] for very short NLSP lifetimes up to mQ= 800 eV/c^. For the 
MSUGRA case the best limit from the fiR has been used. The maximum excluded value 
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of 1^1^= 96.5 GeV/c^ was achieved for 171^= 150 eV/c^. For itLq > 150 eV/c^ the best 
lower mass limit was set by the stable heavy lepton search. Ir masses below 40 GeV/c^ 
were excluded by LEPl data [44]. In the case of Ir degeneracy, this limit improved to 
43 GeV/c^. 

5.2 Neutralino pair production 

Limits for ncutrahno pair production cross-section were derived in the fi NLSP and 
slcptons co-NLSP scenarios for each (m^o,m^-) combination for rriQ below a few eV/c^. 
For the fi NLSP case the combination took into account the results from the LEP runs 
from 1996 (for > 161 GeV) to 2000 [1,32,33]. The hmits for the production cross- 
section allowed some sectors of the (m^o, m^-) space to be excluded. In order to exclude as 
much as possible of the mass plane, the results from two other analyses were taken into 
account. The first is the search for slepton pair production in the context of MSUGRA 
models. From the experimental lower limit on the mass of the fi set by that search, 
rufi > 82.5 GeV/c^ [12], it can be concluded that at least neutrahno masses below the 
lower bound for m-f^ plus the tau mass are excluded. The second is the search for lightest 
neutralino pair production in the region of the mass space where Xi is the NLSP [13] 
(the region above the diagonal line in Figure 10, i.e. ruf-^ > m^o). Within this zone, the 
neutralino decays into a gravitino and a photon. 

As an illustration. Figure 10 presents the 95% CL excluded areas for uIq < 1 eV/c^ in 
the m^o vs nif^ plane for the fi NLSP scenario and for different values of the num- 
ber of messenger generations (n). The positive-slope dashed area is excluded by the 
neutralino pair production search. The negative-slope dashed area is excluded by the 
analysis searching for neutralino pair production followed by the decay Xi ~^ Gj. The 
point-hatched area is excluded by the direct search for slepton pair production within 
MSUGRA scenarios. The neutralino mass range is plotted up to almost the kinematic 
limit. 

5.3 Chargino pair production 

Limits on the production cross-section for chargino pairs were derived for each 
{mQ,mi,m^+) combination for uIq below a few eV/c^. Figure 11 (a) shows, as an 
example, the 95% CL upper limit on the chargino pair production cross-section at 
^/s = 208 GeV as a function of and mj^ after combining with the results of the 
searches for large impact parameter and kink at lower energies using the likelihood ratio 
method [43], for ttIq = 100 eV/c^. The limits on the chargino pair production cross- 
section were used to exclude areas within the {Tn^+ , m^) plane for different domains of 
the gravitino mass, combining results from all the centre-of-mass energies from 183 GeV 
to 208 GeV [1,32,35]. Figure 11 also shows the regions excluded at 95% CL in the 
(m^+,mfj plane, for the fi NLSP scenario (b), and in the {m-+,mj^) plane, for the 

Ir co-NLSP scenario (c). The positive-slope area is excluded for itlq > 1 eV/c^. This 
limit is also valid for smaller masses of the gravitino because they lead to the same fi- 
nal state topologies. The negative-slope area is only excluded for ttIq > 100 eV/c^. The 
areas below mf^= 82.5 GeV/c^ in the fi NLSP scenario, and below ruf^— 88 GeV/c^ in 

the Ir co-NLSP scenario, are excluded by the direct search for slepton pair production 
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in MSUGRA models [12]. These lower bounds on and mj^ from MSUGRA mean a 
limit on the m-± at 82.5 GeV/c^ and 88 GeV/c^, respectively. In Figure 11 the exclusion 
regions given by the search for charginos in GMSB scenarios start at those values. The 
area of Am < 0.3 GeV/c^ is not excluded because in this region the charginos do not 
decay mainly to fi and Ur, but to W and G. Thus, if Am > 0.3 GcV/c^, the chargino mass 
limits are 100 GeV/c^ for aU m^, and 102 GeV/c^ for m^ > 100 eV/c^, in the n NLSP 
scenario. In the sleptons co-NLSP scenario the limits are 96 GeV/c^ for rriQ = 1 eV/c^, 
and 102 GeV/c^ for > 100 eV/c^. The chargino mass limit decreases with decreasing 
rrif^ because in scenarios with gravitino LSP small stau masses correspond to small sneu- 
trino masses (both are proportional to A) and hence to smaller production cross-sections 
due to the destructive interference between the s- and t-channels. 

It should be noticed that within the parameter space studied here, the lightest chargino 
is at least 40% heavier than the lightest neutralino. Thus, for small gravitino masses 
the search for neutralinos implies a lower limit on the lightest chargino of 130 GeV/c^. 
Neutralinos are not directly searched for in heavier gravitino mass regions and therefore 
for this range of gravitino masses the hmit of 102 GeV/c^ for the chargino mass remains 
valid. 

5.4 Heavy stable charged particle pair production 

The results presented in section 4.4 were combined with previous DELPHI results 
in this channel [1,36,37], and cross-section limits were derived within MSSM models as 
indicated in Figure 12. From the intersection points with the predicted cross-sections for 
smuon or staus in the MSSM, left (right) handed smuons and staus can be excluded up 
to masses of 97.7(97.4) GeV/c^ at 95%CL. No limits arc given on selectrons here because 
the cross-section can be highly suppressed by an additional t-channel sneutrino exchange 
contribution. 

5.5 Sgoldstino production 

No excess of events nor clear evidence of an anomalous production of events with 
monochromatic photons is observed in either of the two channels. Therefore a limit on 
the cross-section of the new physics reaction contributing to the two topologies was set. 

The number of detected events, the background rate and the detection efficiency de- 
pend on the (f) mass. In addition, when the expected total width for the same m<^ value 
is comparable or larger than the experimental resolution, the data were compared with 
the background events in a region corresponding to 80% of the signal area. As a conse- 
quence, the limit on the signal cross-section depends on m^ and \/F. Furthermore, to 
take into account the different sensitivities of the two analysed channels, the hkelihood 
ratio method was used [43]. Since the expected branching ratio and total width depend 
on the mass parameters as explained above, the 95% CL cross-section limit was computed 
as a function of m<^ and Vf for the two sets of parameters listed in Table 1, and is shown 
in Figure 13. By comparing the experimental limits with the expected production cross- 
section, it is possible to determine a 95% CL excluded region on the parameter space as 
shown in Figure 14. As explained in [18], to keep the particle interpretation the total 
width r must be much smaller than m^ and therefore the region with F > 0.5m^ was 
not considered. In Figures 13 and 14 the sgoldstino mass range starts at 10 GeV/c^ to 
avoid all theoretical subtleties connected with the non-perturbative aspects of the strong 
interaction. The upper value of the mass range approaches the kinematic Al limit. 
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5.6 Limits on the GMSB parameter space 

Finally, all these results can be combined to produce exclusion plots within the 
(tan /3, A) space. The corresponding parameter space was scanned as follows: 1 < n < 4, 

5 TeV < A < 90 TeV, 1.1 < M/A < 10^ 1.1 < tan/3 < 50, and sign{^) = ±1. As an 
example. Figure 15 shows the zones excluded for n =1 to 4 for ttIq < 1 eV/c^, which cor- 
responds to the NLSP decaying at the main vertex. The shaded excluded. The 
areas below the dashed lines contain points of the GMSB parameter space with Xi NLSP. 
The areas to the right (above for n = 1) of the dashed-dotted lines contain points of the 
GMSB parameter space where sleptons are the NLSP. It can be seen that the region of 
slepton NLSP increases with n. The contrary occurs to the region of ncutralino NLSP. 

The shaded areas below the dashed lines arc excluded by the search for neutralino 
pair production followed by the decay Xi G'j (acoplanar photons search). However, 
the sohd and almost vertical hne that sets the lowest hmit to the parameter A, is defined 
by the combination of the search for acoplanar photons, right handed selectrons with 
sleptons NLSP and ncutralino pair production with stau NLSP. Finally, the search for 
stau NLSP improves the exclusion limit represented by the solid line with positive slope. 

Theoretical arguments constrain the value of tan/3 to be between ~ 1.2 and ~ 65 [8]. 
A lower hmit is set for the variable A at 17.5 TeV. 

6 Summary 

Lightest neutralino, slepton and chargino pair production were searched for in the 
context of light gravitino models. Two possibilities were explored: the fi NLSP and the 
sleptons co-NLSP scenarios. No evidence for signal production was found. Hence, the 
DELPHI collaboration sets lower limits at 95% CL for the mass of the function 
of the slepton masses, and lower mass limits for the sleptons in all the gravitino mass 
range. The limit on the chargino mass is 100 GeV/c^ for all ttiq in the fi NLSP scenario 
and 96 GeV/c^ in the sleptons co-NLSP scenario. All these results were combined to set 
limits on the GMSB parameter space. Combining all the data up to 208 GcV, a lower 
limit is set for the variable A at 17.5 TeV. Similar results have been obtained by the 
ALEPH experiment concerning all these searches [45]. 

Mass hmits for heavy stable charged particles were also derived within the MSSM. For 
these particles the DELPHI collaboration sets lower mass limits at 95% CL for the left 
(right) handed stau and smuon at 97.7 (97.4) GeV/c^. 

Finally, cross-section and mass limits were derived for sgoldstinos at 95% CL since no 
evidence of an anomalous production of events with monochromatic photons was observed 
in either of the two channels. 

Acknowledgements 

We are greatly indebted to our technical collaborators, to the members of the CERN- 
SL Division for the excellent performance of the LEP coUider, and to the funding agencies 
for their support in building and operating the DELPHI detector. 
We acknowledge in particular the support of 

Austrian Federal Ministry of Education, Science and Culture, GZ 616.364/2-III/2a/98, 
FNRS-FWO, Flanders Institute to encourage scientific and technological research in the 



21 



industry (IWT), Belgium, 

FINEP, CNPq, CAPES, FUJB and FAPERJ, Brazil, 
Czech Ministry of Industry and Trade, GA CR 202/99/1362, 
Commission of the European Communities (DG XII), 
Direction des Sciences de la Matierc, CEA, France, 

Bundesministerium fiir Bildung, Wissenschaft, Forschung und Technologic, Germany, 
General Secretariat for Research and Technology, Greece, 

National Science Foundation (NWO) and Foundation for Research on Matter (FOM), 
The Netherlands, 
Norwegian Research Council, 

State Committee for Scientific Research, Poland, SPUB-M/CERN/PO3/DZ296/2000, 

SPUB-M/CERN/PO3/DZ297/2000, 2P03B 104 19 and 2P03B 69 23(2002-2004) 

JNICT- Junta Nacional de Investigagao Cienti'fica e Tecnologica, Portugal, 

Vedecka grantova agentura MS SR, Slovakia, Nr. 95/5195/134, 

Ministry of Science and Technology of the Republic of Slovenia, 

CICYT, Spain, AEN99-0950 and AEN99-0761, 

The Swedish Natural Science Research Council, 

Particle Physics and Astronomy Research Council, UK, 

Department of Energy, USA, DE-FG02-01ER41155. 



22 



References 

[1] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B503 (2001) 34. 
[2] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B494 (2000) 203. 
[3] S. Dimopoulos, S. Thomas, J.D. Wells, Nucl. Phys. B488 (1997) 39. 
[4] S. Ambrosanio, CD. Kribs, S.P. Martin, Phys. Rev. D56 (1997) 1761. 
[5] G.F. Giudice, R. Rattazzi, Phys. Rep. 322 (1999) 419 and Phys. Rep. 322 (1999) 
501. 

[6] D.A. Dicus and S. Nandi, Phys. Rev. D56 (1997) 4166. 
[7] H. Pagels and J.R. Primack, Phys. Rev. Lett. 48 (1982) 223; 
S. Weinberg, Phys. Rev. Lett. 48 (1982) 1303; 

T. Asaka, K. Hamaguchi and K. Suzuki, Phys. Lett. B490 (2000) 136. 
[8] S. P. Martin, hep-ph/9709356. Perspectives on Supersymmetry, ed. by G. Kane, 

Worldscientific, Singapore (1998) pp 1-98. 
[9] S. Dimopoulos, M. Dine, S. Raby, S. Thomas and J. D. Wells, Nucl. Phys. Proc. 

Suppl. A52 (1997) 38. 
[10] E. Calzetta, A. Kandus, F. D. MazziteUi and C. E. M. Wagner, Phys. Lett. B472 

(2000) 287. 

[11] A. Bartl et. al., Z. Phys. C73 (1997) 469. 

[12] EP Paper 302/Draft 2 ( 9 July 2002) (To be submitted to Eur. Phys. J. C). 
[13] DELPHI Collaboration, P. Abreu et al, Eur. Phys. J. C17 (2000) 53; 
DELPHI Collaboration, P. Abreu et a/., Eur. Phys. J. C6 (1999) 371. 

[14] B. Abbott et al, Phys. Rev. Lett. 80 (1998) 442. 

[15] K. Cheung, D. A. Dicus, B. Dutta, S. Nandi, Phys. Rev. D58 (1998) 15008. 

[16] H. Dreiner, hep-ph/9707435. Perspectives on Supersymmetry, ed. by G. Kane, 

Worldscientific, Singapore (1998) pp 462-479. 
[17] DELPHI Collaboration, P. Abreu et al, Phys. Let. B485 (2000) 95. 
[18] E. Perazzi, G. Ridolfi and F. Zwirner, Nucl. Phys. B574 (2000) 3. 
[19] DELPHI Collaboration, P. Aarnio et al, Nucl. Instr. and Meth. 303 (1991) 233. 
[20] DELPHI Collaboration, P. Abreu et al, Nucl. Instr. and Meth. 378 (1996) 57. 
[21] The DELPHI Trigger Group, A. Augustinus et al. The DELPHI Trigger System at 

LEP2 (to be submitted to Nucl. Instr. and. Meth. A). 
[22] T. Sjostrand, Comp. Phys. Comm. 39 (1986) 347; 

T. Sjostrand, PYTHIA 5.6 and JETSET 7.3, CERN-TH/6488-92. 
[23] DELPHI Collaboration, P. Abreu et al, Z. Phys. C73 (1996) 11. 
[24] SUSYGEN 2.20, S. Katsanevas and S. Melachroinos in Physics at LEP2, 

CERN 96-01, Vol. 2, p. 328 and http://lyoinfo.in2p3.fr/susygen/susygen.html; 

S. Katsanevas and P. Moravitz, Comp. Phys. Comm. 122 (1998) 227. 
[25] F.A. Berends and R. Kleis, Nucl. Phys. B186 (1981) 22. 
[26] S. Jadach, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 79 (1994) 503. 
[27] S. Jadach, W. Placzek, B.F.L. Ward, Phys. Lett. B390 (1997) 298. 
[28] F.A. Berends, R. Pittau, R. Kleiss, Comp. Phys. Comm. 85 (1995) 437. 
[29] J. Fujimoto et al, Comp. Phys. Comm. 100 (1997) 128. 
[30] T. Alderweireld et al, in LEP2 MCWorkshop, CERN-2000-009, p 219. 
[31] F.A. Berends, P.H. Davcrveldt, R. Kleiss, Comp. Phys. Comm. 40 (1986) 271, 

Comp. Phys. Comm. 40 (1986) 285, Comp. Phys. Comm. 40 (1986) 309. 
[32] DELPHI Collaboration, P. Abreu et al, Eur. Phys. J. C16 (2000) 211. 
[33] DELPHI Collaboration, P. Abreu et al, Eur. Phys. J. C7 (1999) 595. 
[34] DELPHI Collaboration, P. Abreu et al, Eur. Phys. J. C6 (1999) 385. 



23 



[35] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B466 (1999) 61. 
[36] DELPHI Collaboration, P. Abreu et al., Phys. Lett. B478 (2000) 65. 
[37] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B444 (1998) 491. 
[38] S. Catani, Phys. Lett. B269 (1991) 432. 

[39] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B479 (2000) 129. 
[40] DELPHI Collaboration, P. Abreu et al, Eur. Phys. J. C19 (2001) 29. 
[41] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B462 (1999) 425. 
[42] D. A. Dicus, B. Dutta, S. Nandi, Phys. Rev. D56 (1997) 5748; 

D. A. Dicus, B. Dutta, S. Nandi, Phys. Rev. Lett. 78 (1997) 3055. 
[43] A.L. Read, Modified Frequentist Analysis of Search Results (The CLs Method), 

CERN-OPEN-2000-205, 81. 
[44] Particle Data Group, D.E. Groom et al, Eur. Phys. J. C15 (2000) 1. 
[45] ALEPH Collaboration, A. Heister et al, CERN EP 2002-021, hep-ex 0203024, (subm. 

to Eur. Phys. J.). 



24 
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Figure 1: r mean decay length (L = ctjP) as a function of the gravitino mass for 
different f mass values. 




Figure 2: Sketch illustrating the reconstruction of a secondary vertex in the plane per- 
pendicular to the beam direction. All the radii are measured with respect to the beam 
spot (BS). 
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Figure 3: (a) Angle between the directions defined by the hadronic vertex and the recon- 
structed vertex w.r.t. beam spot, (b) angle between the tracks of the kink, and (c) angle 
between the electromagnetic shower and the missing momentum. Dots are real data, 
cross-hatched histogram is the SM background and blank histogram is the simulated 
signal 



60 GeV/c2, L = 50 cm and = 208 GeV). 
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Figure 4: Efficiency as a function of the decay radius for a sample of staus with L = 50 
cm and = 208 GeV. 
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Figure 5: Acollinearity distribution for real data minus expected SM background 
(cross-hatched histogram), compared to cosmic ray events shown in dark grey. The 
cuts on this variable are shown with arrows. 
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Figure 6: distribution for data (dots), simulated SM backgrounds (grey his- 

togram) and simulated signal of rrif^ = 90 GeV/c^ and itIq = 25 eV/c^ at y/s = 206 GeV 
and a boosted mean decay length ~ 1 cm (white histogram, in arbitrary scale) after all 
other cuts applied by the small impact parameter search. 
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Figure 7: (a) Normalised energy loss as a function of the momentum after the gas veto 
for the 208 GeV data, (b) Measured Cherenkov angle in the liquid radiator as a function 
of the momentum after the gas veto: if four photons or less were observed in the liquid 
radiator, the Cherenkov angle was set equal to zero. The areas labelled (I), (II) and 
(III) indicate the selection criteria explained in the text. Open circles are data. The 
small filled circles indicate the expectation for a 95 GeV/c^ mass signal with charge ±e, 
resulting in a large dE/dx (upper plot) and no photons (except for a few accidental rings) 
in the liquid Cherenkov counter (lower plot). The solid lines with a mass signal value 
indicate the expectation for heavy stable sleptons. 
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Figure 8: a) Photon recoil mass spectrum for 777 candidates (points) and the expected 
background (histogram). The average number of entries per event in the data is 2.3. 
The bin size takes into account the experimental mass resolution and the expected signal 
width, b) Photon recoil mass spectrum for '~^gg candidates (points) and the expected 
background (histogram) . 
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Figure 9: Exclusion regions in the {mQ,mfJ (a) and {mQ.mi^ (b) planes at 95% CL 
for the present analyses combined with the Stable Heavy Lepton search and the search 
for / in gravity mediated models (MSUGRA), using all DELPHI data from 130 GeV 
to 208 GeV centre-of-mass energies. The dashed line shows the expected limits for the 
impact parameter and kink searches. 
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Figure 10: Exclusion regions in the {m^o^rrif^) plane at 95% CL for the neutralino analysis 
(positive-slope dashed area) combined with the search for neutralino pair production 
followed by the decay Xi G'-y (negative-slope dashed area), and the direct search for 
fi pair production within MSUGRA scenarios (point-hatched area), using all DELPHI 
data from 161 GeV to 208 GeV centre-of-mass energies. 
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Figure 11: (a) Limits in picobarn on the lightest chargino pair production cross-section 
at 95% CL. Limits are shown as a function of mj and m-+ for ttIq = 100 eV/c^. Areas 
excluded at 95% CL in the (m-+,mfj plane (b) and {m^+,mi^ plane (c). The posi- 
tive-slope hatched area is excluded for all gravitino masses. The negative-slope hatched 
area is only excluded for itlq > 100 eV/c^. The grey area is excluded by the search for 
stau pair production in gravity mediated supersymmetry breaking models. Both plots 
have been obtained using all DELPHI data from 183 GeV to 208 GeV centre-of-mass 
energies. 
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Figure 12: Predicted production cross-section for left and right handed stable smuons 
(staus) as a function of the particle mass. The cross-section limit indicated in the figure 
has been derived using all DELPHI data between 130 and 208 GeV. 
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Figure 13: Cross section upper limit (pb scale on the right) at the 95% CL as a function 
of and ^/F for the two sets of parameters of Table 1 and using all DELPHI data 
between 189 and 208 GeV. 
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Figure 14: Exclusion region at the 95% CL in the (m^,V^) plane using all DELPHI 
data between 189 and 208 GeV. Regions 1 and 2 are the excluded regions once the set of 
parameters (labelled with 1 and 2 respectively) of table 1 are applied. 
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Figure 15: Shaded areas in the (tan /3, A) plane are excluded at 95% CL. The areas below 
the dashed lines contain points of the GMSB parameter space with Xi NLSP. The areas 
to the right (above for n = 1) of the dashed-dotted lines contain points of the GMSB 
parameter space were sleptons are the NLSP. 



